Mapping small inland wetlands in the South-Kivu province by integrating optical and SAR data with statistical models for accurate distribution assessment

There are several techniques for mapping wetlands. In this study, we examined four statistical models to assess the potential distribution of wetlands in the South-Kivu province by combining optical and SAR images. The approach involved integrating topographic, hydrological, and vegetation indices into the four most used classifiers, namely Artificial Neural Network (ANN), Random Forest (RF), Boosted Regression Tree (BRT), and Maximum Entropy (MaxEnt). A wetland distribution map was generated and classified into 'wetland' and 'non-wetland.' The results showed variations in predictions among the different models. RF exhibited the most accurate predictions, achieving an overall classification accuracy of 95.67% and AUC and TSS values of 82.4%. Integrating SAR data improved accuracy and precision, particularly for mapping small inland wetlands. Our estimations indicate that wetlands cover approximately 13.5% (898,690 ha) of the entire province. BRT estimated wetland areas to be ~ 16% (1,106,080 ha), while ANN estimated ~ 14% (967,820 ha), MaxEnt ~ 15% (1,036,950 ha), and RF approximately ~ 10% (691,300 ha). The distribution of these areas varied across different territories, with higher values observed in Mwenga, Shabunda, and Fizi. Many of these areas are permanently flooded, while others experience seasonal inundation. Through digitization, the delineation process revealed variations in wetland areas, ranging from tens to thousands of hectares. The geographical distribution of wetlands generated in this study will serve as an essential reference for future investigations and pave the way for further research on characterizing and categorizing these areas.


Abbreviations
Artificial Neural Network (ANN), Logistic Regression Model (LRM), and Maximum Entropy (Maxent) have been utilized, with Random Forest (RF) and SVM showing promising results in terms of accuracy.These models have also been tested for wetland distribution and other applications, proving effective 35,36 .Overall, both DT, SVM, ANN, RF, BRT, KNN classifiers are the most famous worldwide classification algorithms used for wetland mapping 37 .Despite abundant data and tools, there still needs to be more knowledge regarding which models can accurately assess the practical distribution map of small inland wetlands and delineate them.Debates persist regarding the choice of models and the types of images or indices to employ for this purpose.However, advancements in technology have significantly contributed to the field of Geographic Information Systems (GIS) and Remote Sensing (RS), enabling the modeling and prediction of wetland ecosystems at both small and macro scales, as well as the assessment of distribution factors.Modern geo-statistics and techniques integrated into GIS tools facilitate efficient modeling of wetland distribution, while RS provides valuable imagery for detecting, digitizing, and estimating wetland distribution.These advancements have enhanced our understanding and capability to assess wetland distribution accurately, contributing to effective wetland management and conservation efforts.
To address the need for identifying and delineating these low-lying wetlands in the eastern region of the DRC, we propose an approach that combines optical images from the Sentinel satellites with Synthetic Aperture Radar (SAR) images.The methodology draws upon previous research by Mwita et al. 30,38,39 , and Garba et al. 32 , with modifications such as replacing Landsat images with Sentinel (1 and 2) and ALOSPALSAR.Additionally, we evaluate the performance of four widely used statistical classifiers in wetland mapping.
These statistical models have been extensively studied and proven to enhance the accuracy of wetland distribution predictions 40,41 .Over the past decade, these techniques have garnered significant attention in ecosystem modeling and forecasting thanks to their ability to improve predictive capabilities.One key advantage of these mathematical models is their utilization of different types of independent and dependent variables, including categorical and quantitative variables.This versatility extends their applicability beyond wetland mapping and serves as a valuable reference for researchers in various fields of science, such as biology, sociology, and agronomy.

Study objectives
The overall objective of this study is to contribute to the identification and study of wetlands in the Democratic Republic of the Congo (DRC) by developing improved methods for mapping small inland wetlands.Specifically, this study aims to achieve the following objectives: (1) Identify the critical explanatory variables derived from remote sensing data, including Sentinel-1 and Sentinel-2, as well as ALSOPALSAR, and field data, that are relevant for modeling the distribution of small wetlands in eastern DRC; (2) Evaluate the capabilities of singledate Sentinel optical data and their combination with Synthetic Aperture Radar (SAR) data for mapping small wetlands in the South-Kivu province.This evaluation will include assessing the accuracy of these mapping methods and identifying potential errors; (3) Digitize and characterize the identified wetlands, including analyzing their morphological characteristics such as area and perimeter; (4) Discuss the strengths and limitations of the mapping methods employed in this study, providing an overview of the advantages and challenges associated with integrating optical, topographic, and SAR indices, as well as using novel classifiers, to achieve accurate mapping results.
To address these objectives, it is hypothesized that integrating optical, topographic, and SAR indices with new classifiers will result in a more accurate method for mapping small wetlands.Additionally, it is suggested that fully polarimetric SAR imagery can provide valuable information about surface scattering mechanisms, allowing for a more precise distinction of small wetlands.Furthermore, including SAR data and novel vegetation indices is expected to improve the mapping process.Finally, it is anticipated that the delineation of wetlands after digitization will reveal various areas with distinct morphological characteristics.By addressing these research hypotheses, this study will enhance our understanding of wetland mapping methodologies in the context of the DRC, with a specific focus on small wetlands.

An overview of the study area
The South-Kivu province is situated in eastern DRC and shares borders with Rwanda, Burundi, and Tanzania.Covering an area of approximately 64,791 km 2 , it is one of the 26 provinces of the Democratic Republic of the Congo.South-Kivu accounts for around 2.73% of the total land area of the country 42,43 .The province is home to approximately 6.2 million people, with 47% of the population residing in rural areas 44 .Geographically, the province is located between latitudes 1.5836° and 5.0103° South and longitudes 26.8106° and 29.3890°East (Fig. 1).
The province is divided into eight territories, namely (Shabunda, Kalehe, Uvira, Walungu, Kabare, Idjwi, Mwenga, and Fizi) including the city of Bukavu (Fig. 1), and each territory in municipalities referred to as "Groupement".The province has a variety of landscapes ranging from low to very high altitudes.The elevation fluctuates from ~ 512 to 3464 m above the sea level (m.a.s.l) and decreases from east to west.The Congolese "Cuvette Centrale" begins in the Shabunda and Mwenga territories, whereas reliefs and valleys, including mountains and the Mitumba Chain Mountains, characterize the eastern and central territories.On the left, the Ruzizi plain is a broad plain extending into the Fizi territory, in the territory of Walungu and Uvira, as well as the highland favorable to the growth of a diverse type of wetlands and especially swamps, marshes, and peats 44 .Various wetlands occur due to these physical conditions (Figs. 2 and 4).
South-Kivu province is characterized by a humid tropical and equatorial climate.According to Balasha et al. 6 , the mean annual rainfall of ~ 1500-1800 mm each year and an average annual temperature varies between 11° and 25 °C (Supplementary material S1) are observed in the province.Seven primary soil categories predominate in the South-Kivu mainly Haplic Acrisols, Humic Cambisols, Humic Ferrasol, Luvic Phaeozems, Mollic  www.nature.com/scientificreports/Fluvisols, Vertisols, and Gleyic Solonchaks 45,46 .Histisols (organic soil) are observed in some territories in small areas.The peat is used as energy fire material for food cooking in small quantities using small bricks.These bricks can replace a sizable quantity of charcoal or wood.In the South Kivu, peat is hand-harvested in the Kakonda peat (in Kabare), Hogola (at Nyangezi) in the Chiherano, and Kachandja peats (in Walungu).The hydrography is extensive and thick, with several small sources, large rivers (such as the Ruzizi, Elila, Ulindi, Itombwe, and Lwama), as well as lakes.Additionally, the province is a part of the "African Great Lakes" region.It has Lake Kivu, Tanganyika, and a few more small lakes and ponds (such as Lubwe in Itombwe).The vegetation comprises highland forests, herbaceous savanna, wooded swamps, and dense forest.The land cover map was taken from the 2018 version of the European Space Agency's (ESA) Climate Change Initiative Land Cover Project (CCI-LC) (cds.climate.copernicus.eu)product.The CCI-LC identified 21 land cover classes in South-Kivu (Supplementary material S2), which can initially be categorized using the FAO's land cover classification system (LCCS).The province has protected areas, mainly the Kahuzi-Biega National Park (KBNP), the Natural Reserve of Itombwe (NRI), Maniema, South-Masisi, the Mont Kabobo (National Park of Ngamikka), and the Luama-Kivu Hunting Reserve 47 .

Materials and methods
A complete overview of the methodological flowchart can be found in Fig. 3.The method used is discussed in detail in the following subsections comprising data pre-processing, classification, training data, and statistical models used for prediction and validation.In this study, we adapted the methodology suggested by Mwita et al. 38 for small inland wetlands of East Africa; and revised it according to the one proposed by Adeli et al. 39 , LaRocque et al. 30 , and Garba et al. 32 .We combined Sentinel-1 and -2 and ALOS-2 PALSAR data; four statistical classifiers were used.The first step was downloading satellite images.These images were obtained from the ESA and JAXA websites (https:// www.eorc.jaxa.jp/ ALOS/ en/ index_e.htm).Once obtained, the processing follows as the second step.It consisted of optical image extraction, correction, and band merge.The corrected bands were clipped to the study area (South-Kivu shapefile obtained from RGC: http:// rgc.cd/), and indices were calculated.The indices calculation is presented in Supplementary material S3.The indices formulas were extracted from (https:// custom-scrip ts.senti nel-hub.com/ custom-scrip ts/ senti nel-2/ index db/).In total, 27 indices were used.

Distribution data
The term "wetland" used in this study is based on the definition suggested by Steinbach et al. 48, Chuma et al. 5 , and Amler et al. 7 , and which is currently used in eastern Africa and adapted from the Ramsar definition 49 : 'wetlands www.nature.com/scientificreports/are areas of marsh, swamp, inundated valleys, peatland or water, whether natural or artificial, permanent or temporary, with static or flowing water, the depth of which at low tide does not exceed 6 m' .Conceptually, the "potential distribution, existing wetlands" method has served as the foundation for this study.It concentrated on prospective wetlands, representing wetlands' size at the most significant level before human interference.In other words, regions where water-related ecosystems are most likely to develop, are highlighted through prospective wetland mapping.Archive points and polygon files constitute the presence data, which were collected by scientists, environmental local non-government organizations (NGOs), and mostly during our field works from September 2020 to August 2022.Other visual fieldworks of well-known sites were made using drone images and delineated in ArcGIS and Google Earth following the digitalization process.In total, 550 samples shapefiles were considered (Supplementary data S4).For inaccessible wetlands, the location's geographic coordinates were taken from the edge and placed in the center once the image was obtained.More points were taken into consideration for the wetlands that had a large area.The total sample (shapefiles) was split into two datasets; 70% of plots (385 shapefiles) were used as training samples, and the remaining (30%: 165) as validation samples.Those points were considered "presence samples".The "Absence" data comprised 2300 shapefiles extracted from the RGC comprising cities and main towns, schools, hospitals, villages, airports, farms, tree plantations, woodland, etc.These ground-truth surveys were conducted during the same period as the image acquisition.

Data on environmental variables
Since wetlands are characterized by specific topography, vegetation, and hydrology, topographic, vegetation, and hydrological indices were used during the mapping and delineation process.Topographic variables were derived from the digital elevation model (DEM) from ALOS-PALSAR (12.5 m resolution).This comprises the elevation (m), the slope (%), land aspect, curvature, Topographic Witness index (TWI), also known as the compound topographic index (CTI), estimated following the formula: TWI = (αtanb) where α is the local upslope area draining through a certain point per unit contour length and b: the local slope in radians.TWI is related to soil moisture influencing rapid runoff and flash floods 50 .It was calculated in ArcGIS 10.7 Esrti-TM.The ALOSPALSAR DEM was first projected, and the flow direction and flow accumulation were calculated.TWI was calculated as the ratio between the Flow accumulation and the slope.For Sentinel-1, the horizontal transmit and horizontal receive polarization (HH), horizontal transmit and vertical polarization (HV), and vertical transmit and vertical receive polarization (VV); the ratio γ =|HH| 2 |VH| 2 and η =|VV| 2 |VH| 2 were also calculated and integrated.For ALOS-2 PALSAR data, three features were integrated: HH, HV, and HH/HV (ρ =|HH| 2 |VV| 2 ).They were corrected using Freeman-Durden (following surface scattering and double-bounce scattering), Cloude-Pottier (polarimetric decompositions and the compact polarimetric simulations) comprising alpha, beta, gamma, and lambda, multipolarizations, dual polarizations, and polarimetric decompositions 51 .The SAR data preprocessing comprised speckle reduction, terrain correction, and geocoding following steps developed by Veci 52 , Foumelis et al. 53 , and Braun 54 .Since the province typically has two seasons, one dry season lasting from May to August and another dry season for the remaining months, two images were used, comprised of those taken in July (at the height of the dry season) and one in November (during the high rainy period).
This study considered hydrogeomorphic (HGM) wetland types: riverine, depression, slope, flat, and lacustrine fringe.Wetlands were mapped by modeling groundwater using environmental variables.The standard algorithms implemented steps in S1 Toolbox software include applying orbit files, removing low-intensity noise and invalid data on scene edges, removing thermal noise, radiometric calibration, and orthorectification 51,54 .The correction was made in the Sentinel Application Platform (SNAP) for optical data, starting with Rayleigh correction (computing the bottom of Rayleigh reflectance bands).The image resolution (in m) was resampled at 10 m.The sea level pressure and ozone were maintained at 1013.25 (in hPa) and 300 (in DU).Finally, radiance-to-reflectance conversion was used before indices calculation 55 .The dataset of exploratory variables included 27 environmental variables divided into 7 topographic and Hydrogeomorphology (HGM) variables, 8 SAR, and 12 vegetation indices.The description, formula, and source of these indices are presented in the Supplementary material S3.The radar backscattering was made following Sigma-nought (σo).It is referred to as the radar backscatter per unit area (m 2 /m 2 ), expressed in decibels (dB).The standard formula used to calculate σo: σo = 10 * log10(DN2) + K, where DN is the image pixel digital number measured in the SAR amplitude image, and K is a calibration factor that varies depending on the SAR sensor and processor system used.As we use both ALOSPALSAR and Sentinel-1, the factor was fixed to − 83.0 dB and 0 dB, respectively.
Other ancillary information was used, including GPS points, on-site images obtained from drone missions (Fig. 4), field notes on dominant vegetation, accessibility from roads, villages, or markets, and land use.The GPS points were inserted into ArcMap and google earth, and then the boundary delineation was conducted using high very resolution images.The Universal Transverse Mercator (UTM), zone 35S, was used for projection, while the Geographic Coordinate System (latitude-longitude) WGS1984 was maintained for points.Three non-wetland classes: deep water, urban, and upland, were obtained and merged into one class called 'no-wetland' .

Model selection and construction
Four spatial statistic models were used comprising Artificial Neural Network "ANN" 21 , Boosted Regression Tree "BRT" 56 , Random Forest "RF" 11,56 and Maximum Entropy 'Maxent' 35 .It is a computing system inspired by the biological neural networks that make up animal brains, known as neural networks (NNs) or neural nets.
(a) ANN is a widely used machine learning (ML) algorithm that can work for big data analysis.The multilayer feed-forward feature is the primary form of neural networks.ANN includes several neurons or nodes that function in parallel to convert the input data into output types.Typically, ANN consists of three-layer types, namely (i) the input, (ii) the hidden, and (iii) the output layers.Depending on the specific application in a network, each layer has some neurons.Each neuron is connected to other neurons in the next consecutive layer by direct links.These links have a weight that represents the strength of an outgoing signal 57 .
(b) The Random Forest (RF) model has been widely used to map this ecosystem 58,59 .According to Rapinel et al. 11 , due to its ability to consider a large number of variables from many sources and low sensitivity to outliers and over-learning, RF models have been more effective for mapping wetlands than other types of models like support vector machines (SVM), maximum likelihood, or decision trees (DT).
It is a non-parametric supervised ML algorithm inside RStudio and R 4.2.1 60 .Because it can manage the significant difference variables and be used to neutralize noisy data, RF has demonstrated its use for classifications with the enormous volumes of data in satellite images.The input variables for the Random Forests were indices generated from the satellite images and sample points (presence and absence).The number of trees used in each Random Forest classification was set at 250, which has been found to be the ideal amount when accuracy and processing speed are considered.To reduce the tree depth, the minimum node size was set at five.The RF algorithm builds numerous bootstrapped, de-correlated random decision trees to categorize a dataset according to the mode.Indeed, when implementing the RF classifier must specify the number of decision trees and randomly chosen variables for dividing the trees.The RF technique was developed with 250 decision trees utilizing the stratified random sample points after multiple iterations and fine-tuning.An ensemble of 10 models was generated for each combination.The number of training samples was set to 5000 and maintained 10 RF models.
(c) Boosted Regression Tree (BRT): To maximize prediction accuracy while revealing information about pertinent variables and their interactions, BRTs adaptively construct several; basic regression-tree models and merge them into a multi-tree model.This differs from conventional regression methods, which result in a single prediction model.With the added benefits of boosting, which makes it possible to model nonlinear functions and improves robustness to data concerns like outliers, BRTs combine the flexibility of regression trees to accept all data types inside a model, including missing or non-independent data.BRTs are described in depth and detail by Berhane et al. 59 .A BRT model needs two crucial inputs: learning rate and tree complexity.While the latter identifies the number of nodes inside each tree and sets the number of variable interactions that are fitted, the former establishes each tree's contribution to the final model 60 .As each tree contributes less to the overall model with a slower learning rate, more trees are in a BRT model.However, building the model requires more computing effort and observations the more trees there are.Averaging 550 trees per model, we used a learning rate of 0.005 and a tree complexity of 5.In order to reduce the inherent stochasticity in each model because of the subsampling and bagging that go into the construction of each tree, an ensemble of 10 models was created www.nature.com/scientificreports/for each combination of inputs and then averaged.Boosted regression trees (BRT) combine the strengths of two algorithms: regression trees (models that relate a response to their predictors by recursive binary splits) and boosting (an adaptive method for combining many simple models to give improved predictive performance).The final BRT model can be understood as an additive regression model in which individual terms are simple trees fitted in a forward, stage-wise fashion.
(d) Maximum of Entropy (Maxent): Currently used for species distribution, Maxent assumes all presence locations of the same type (usually species), and there is little variability in the species' niche preference for locations.It assumes that environmental conditions at extant wetland locations represent the fundamental niche for a particular type of wetland.Maxent cannot adjust for differences in sampling effort type of wetland.In our case, the remotely sensed wetland presence data ensured an unbiased sample, while modeling limitations were collectively addressed using extensive wetland presence coverages to perform validation studies.In this study, we have adapted the methodology developed by Rebelo et al. 61 .We set the number of iterations to 5000, allowing the model enough time to converge.All other Maxent settings were left at their default values.Ten models were also maintained.This technique uses environmental data (here variables) for several "background" sites and known presence locations.The raster files contain variables that were extracted.The result, which was displayed as a probability raster, was reclassified into two classes above-mentioned.
Data were first processed in ArcGIS 10.7 Esri-TM and R 4.2.1 to run the models.The R packages "raster", "rgdal", "virtual species", RStoolbox", "xlsx", and "data.table"help to process the images and the training samples.These packages also helped with data stacking, data frame creation, and raster processing to the same extent, resolution, and projection system.All the data were brought to WGS84/UTM zone 35S.To run the models, the package "MASS", "gbm", for BRT, "dismo", "SDMtune", "NMeval", "virtualspecies", "maxnet" for Maxent, "neuralnet" and "caret" for ANN, "randomForest" and "caret" for RF.For the result presentation, "ggplot2" and "gpairs" helped design graphs.The output files are probability images (ranging from 0 to 1) that were split into two classes, "wetland" and "non-wetland" 1, 35 .

Post classification
The post-classification stage of this study aimed to refine the wetland inventory by removing erroneously classified pixels.Indeed, before conversion, post-classification was made.Three spatial analyses were made; we started first with the "majority filter tool" followed by the "boundary clean tool".The wetland shapefiles were smoothed using the boundary between areas by expanding and shrinking.The small clusters were processed using the "region group tool" 62 .The "polynomial approximation with exponential kernel smoothing algorithm" was finally used to smooth wetland shapefiles.

Accuracy assessment
To determine the level of accuracy of the small inland wetland map, it is necessary to validate the predicted wetted landscape using optical and SAR imageries.The second dataset (30%: 165) of wetlands was used for validation.These wetlands were physically surveyed to explore their nature and water seasonal availability to validate the predicted wetland areas.300 sites were collected from Google Earth Pro images and selected from a field survey using GPS.They were used to assess the overall and local accuracy.The accuracy was assessed using the Area under the curve (AUC), Kappa coefficient, correlation, and True Skill Statistics (TSS).First, we started assessing the "Sensitivity" and "Specificity".In fact, due to their complementarity and the fact that they comprise the core of the Receiver Operator Characteristic "ROC" curve, Specificity and sensitivity are two fundamental metrics in classification 63,64 .Sensitivity and Specificity in this study relate to the accuracy of the "wetland" and "non-wetland" classes producer, respectively.Both concepts have their roots in the literature on ecological presence-absence distribution models, such as species and habitat distribution models 65 , which are comparable to the probability of wetland-occurrence models developed.Here, sensitivity is a term used to define the percentage of presences (wetlands) identified; it measures the degree of omission errors.Specificity is the fraction of accurately anticipated absences (i.e., non-wetland) that occurred.The formula (i) and (ii) are used to ass the two parameters and combined in the confusion matrix as presented below.The two parameters are True Positive Rate (TPR) and False Positive Rate (FPR).Sensitivity is generally used to describe the proportion of observed presences (i.e., wetlands) that are correctly predicted as such and is a measure of the level of omission errors.Specificity, on the other hand, is a measure of errors of commission, representing the proportion of correctly predicted absences (i.e., non-wetland).Sensitivity = α α+δ (i) and Specificity = β β+γ (ii), with n = α + β + γ + δ. α (TP) is the true presence, β (TN): true negative, γ (FP): the false positive and δ (FN): the false negative.
AUC: is a comprehensive assessment of the model's overall performance.It provides a conceivable categorization threshold and is interpreted differently; one is the probability that the model classifies a positive example more carefully than a negative one.AUC is a robust and commonly used measure of predictive model performance.The AUC ranges from 0 to 1.An entirely false prediction model will have an AUC of 0, whereas an accurate prediction model will have an AUC approaching 1, 63,64 .The value 1 specifies that the diagnostic test is perfect, 0.5 stands for a worthless test, and 0 signifies the result is entirely wrong.According to Kanti et al. 65 , the AUC value of 0.90-1.00indicates excellent, 0.80-0.90means good, 0.70-0.80means fair, and 0.60-0.70means poor accuracy level.
Correlation coefficient (r 2 ): The correlation coefficient (r 2 ) is a standardized measure of the predictive accuracy of a model.The formula (iii) was used to assess the r 2 .
Y obs is the wetland, and Y sim is the value from the model.Ȳ obs , Ȳ sim are the average values of Y obs and Y sim .As for AUC, r 2 varies from 0 to 1. 0.90-1.00indicates excellent, 0.80-0.90means good, 0.70-0.80means fair, and 0.60-0.70means poor accuracy level 66 .
Kappa statistic: The Kappa coefficient is more valuable than the overall accuracy as it indicates how the classification rate compares to the likelihood of correctly classifying pixels by chance.The formula (iv) was used to assess the Kappa coefficient.

True Statistic Skill (TSS):
The TSS, often called the Hanssen-Kuipers discriminant, compares the number of correctly classified samples to that of a hypothetically perfect classification while removing those correctly classified samples that could be attributed to a chance agreement.Equation (v) below is used to calculate the TSS, which provides a more objective measurement of classification accuracy comparable to the Kappa coefficient often employed in the literature on remote sensing.

Contribution of environmental variables
The study employed a rigorous variable selection process to identify the most influential variables for the modeling procedure.Of all the variables integrated into the models, 27 were integrated into the modeling process.The process used for variable selection comprises the Pearson correlation (r) calculation followed by the distance calculation (D = 1 − Pearson's r).A two-by-two variables comparison, with D < 0.3 (r > 0.7), was used for variable selection.The tree plot that exhibited strong correlations was generated and provided in the additional data.If two are highly correlated (r > 0.7), only one variable was selected to avoid data multi-collinearity (Supplementary material S5).From 27 indices, only 15 were selected.Indeed, DGM, BBLUE, GNDVI, NDWI, and RVI were highly correlated, and only NDWI was selected.WET, GVMI, and NBR were also correlated, and WET was selected.NDVI was selected among MSAVI 2, TVI, SAVI, and OSAVI.Explanatory variables (indices contribution) values vary from one model to another.Two SAR indices were integrated to assess the "permanent" and "seasonal" effects.The results indicate that for all the wetland and non-wetland classes, the ratio ρ and η for both seasons performed better than the other two Sentinel 1 variables (ratio γ, HH, VV, or VH).Table 1 presents the contribution of each explanatory variable for the four models.

Model prediction
To evaluate the performance of the four models, we started with the accuracy assessment results presented in Fig. 5. First, using a model without SAR data.With 10 experiments for each model, statistical analyses were used for comparison.The AUC varies significantly from one model to another (IC 95 [0.45-0.77]).More variability was observed in the BRT, ANN, and MaxEnt than in the RF models.The decreasing order of the average AUC values is 0.97, 0.96, 0.95, and 0.91 for RF, ANN, MaxEnt, and BRT, respectively.Thus, even if all the models presented high AUC, the best accuracies were achieved for the experiment from RF.For the coefficient of correlation, the same trends were observed with 0.82, 0.67, 0.68, and 0.58 for RF, ANN, MaxEnt, and BRT, respectively.RF and MaxEnt had acceptable values.Only RF (0.96) presented a high value for the Kappa coefficient compared to all the other models.The same trend was observed for TSS (0.82).Considering the TSS as the final model assessment, it is clear that RF (0.84) has a higher and acceptable value than others.Other models, such as ANN and MaxEnt, presented high variabilities in the experiments.TSS values vary from acceptable to unacceptable values (0.28-0.81).The RF model consistently presented significant values followed by ANN, MaxEnt, and BRT for all model evaluation parameters.Even though all the models generally showed high values for AUC and TSS, they presented very high variabilities (see the height of the boxplot) except for the RF.Table 1 presents the contribution of explanatory variables.As above-mentioned, from the 27 variables integrated into the models, only 15 were significantly used, and others were omitted to avoid multi-collinearity.11 variables contributed to RF and ANN models.In contrast, 10 and 9 were integrated for BRT and MaxEnt models.Other variables had minimal contribution (with contribution < 1).Five factors have highly contributed: terrain slope, TWI, GOSSAN, MNDWI, ρ, and η.Topographic indices highly come in the first position.For ANN, ~ 78% of contribution comes from topographic indices, mainly slope (52%), elevation in m (11%), TWI (7%), and Curvature (8%); while for BRT, topography contribution reaches 74% with slope (53%), Aspect (8%), elevation (8%), nNDVI (8%) and TWI (5%).For the MaxEnt model, topographic indices reach 92% of contribution and ~ 5% for vegetation.RF combined more parameters including topography: 68% (slope: 33%, TWI: 17%, Aspect: 6%, elevation: 6%), vegetation (NDWI: 7%, NNDVI: 9%, EVI: 3%).

Prediction of distribution of inland wetlands
Spatial analysis was employed to characterize the "wetland" and "non-wetland" classes by predicting their area (in hectares: ha) and perimeter (in each territory.The results indicate an average wetland coverage of ~ 13.7% (898,690 ha) in the South-Kivu province, with variations observed across different models (Fig. 7).The ANN model predicted a wetland coverage of ~ 14% (967,820 ha), while the MaxEnt model estimated it to be around 15% (1,036,950 ha).The BRT model resulted in a higher wetland coverage of ~ 16% (1,106,080 ha).In contrast, the RF model predicted a lower wetland area at the provincial scale, with a coverage of 10% (691,300 ha).Considering the higher accuracy of the RF model, one could infer that the first two models, which slightly underestimate the wetland surface area, are more aligned with the actual scenario compared to the other two models.
The RF model without SAR data predicts that wetland coverage in the South-Kivu province is approximately 10% of the total area.However, when SAR data is integrated, the proportion increases to 13.5%, closer to the merged model's overall mean value.This suggests that roughly 3.5% (241,955 ha) of the province consists of seasonally flooded wetlands.The wetland distribution was then classified and converted to shapefiles, followed by clipping based on each territory.Figure 8b illustrates the calculation of wetland and non-wetland class areas concerning the territory areas.The surface area and perimeter of wetlands vary across different territories.Fizi has the highest wetland surface area at 30% (47,364 ha), followed by Mwenga at 29% (32,398 ha), Shabunda at 21% (52,743 ha), and Uvira at 19% (5981 ha).Conversely, Idjwi has the smallest percentage of wetlands at 9% (252.9ha), followed by Kalehe at 6% (3075.6 ha), Walungu at 13% (23,400 ha), and Kabare at 16% (31,360 ha).Despite their relatively minor share of the overall area, Shabunda's wetlands cover a larger area than other regions (Fig. 9).The wetland maps generated using the four machine-learning models are presented in Fig. 10a,b.
Based on these maps, particularly Fig. 10b, the western territories are characterized by permanent wetlands, while the eastern wetlands are predominantly seasonally flooded wetlands.The figure also reveals that seasonally flooded wetlands primarily characterize Uvira, Fizi, and Kabare.In terms of proportions, areas with permanently flooded wetlands, with a persistent layer of water on the soil's surface throughout the year, are more prevalent (~ 68%) than areas with seasonally flooded wetlands.The study area description mentions that the South-Kivu 0.8 0.9  province has protected areas.Overlaying the wetland map with the natural reserves (NRs) in South-Kivu, it is evident that most of these wetlands lie outside of these NRs or protected areas.Figure 10 highlights areas with larger wetland surfaces or multiple small wetland areas grouped, which can be considered "wetland complexes."These areas are present in each territory but are particularly prominent in Fizi, Mwenga, and Shabunda.Examples include the complex near Lake Lubwe and the Kalungwe, Kibu, and Mwana rivers in the NRI (Fig. 11) and some in the KBNP.In Fizi territory, the Kilombwe River complex and Lubishako, situated in the Kabobo Mont (NR of Ngamikka), and another small wetland complex near the Nemba and Mulambala rivers in northern Fizi can be observed.Overlapping the main river shapefile reveals that these wetlands are located along major rivers and coastal lakes, such as the Tanganyika in the case of Fizi.In Uvira territory, the region surrounding the Ruzizi River, particularly the small Ruzizi delta, rich in swamps and peat, and the entire length of the river, are abundant in wetlands.Other wetland complexes in the Mwenga territory include the meanderings of the Elila, Semuliki, Ulindi, Nezemere, Kibu, and Kalungwe rivers: the Kilombwe and Lubishako rivers, the Luama complex, and the Kandja (Fig. 9).
In the northern part of Fizi territory, an additional wetland complex is formed by the Sundja, Mulambala, and Nemba rivers outlets, located along the shoreline of Lake Tanganyika.Uvira territory predominantly features wetlands along the Kiliba, Mulongwe, Sange, and Luberizi rivers.These wetlands have undergone significant transformations and are now mainly utilized as inundated valleys for rice cultivation, and in the case of Kiliba, for sugar cane farming.This entire region is commonly referred to as the Ruzizi valley plain.In Walungu territory, wetlands can be found in Kamanyola (along the Ruzizi River), Kaziba chiefdom, Nyangezi, and Walungu-Ciherano axis.The southern part of the Kabare territory consists of small marshlands within enclosed valleys.The largest wetland complexes in the Kabare territory are still within the Kahuzi-Biega National Park (KBNP) and between the territories of Shabunda and Kalehe.Speaking of Shabunda territory, there are extensive wetland complexes along rivers, including notable ones along the Lugulu, Duma, Ulindi, and Kasema Rivers and along the Mosala River.The next step involved using high-resolution images to delineate and characterize many wetlands.Each wetland's total surface area and perimeter were determined by calculating the shapefiles obtained from the digitization process and the conversion of the prediction map.The size of wetlands varied greatly, ranging from a few tens to thousands of hectares.In Ijdwi territory, the wetland ranged from 1.3 to 12.3 ha (IC 95 = [1.4,5.6], IC 95 : confidence interval of 95%).Kalehe territory exhibited wetlands ranging from 2.4 to 620 ha (IC 95 = [4.2,87.8]), while Walungu territory had wetlands varying from 2.5 to 860 ha (IC 95 = [3.7,175.2]).Smaller wetlands primarily characterized these territories.The other territories feature wetlands with areas ranging from 5 to 9860 ha and perimeters ranging from 0.8 to 122.6 km.Large wetlands were observed in Fizi, ranging from 1.5 to 8660 ha with an average of 559 ha, IC 95 = [254, 2600]), and in Uvira, ranging from 1.8 to 2370 ha with an average of 163 ha, IC 95 = [122.6, 378]).For Shabunda and Mwenga, the average wetland sizes were 5 to 2450 ha (IC 95 = [58.8,245.8]  with an average of 85.3 ha, and 6.3 to 1781 ha (IC 95 = [34.3,132.4] with an average of 71.3 ha respectively.In Kabare, wetland areas varied from 4.1 to 3270 ha with an average of 49 ha (IC 95 = [120.7,225.8]) (Fig. 12).Overall, at the provincial scale, the average wetland area and perimeter were ~ 163 ha and ~ 123.6 km, respectively.

Mapping inland wetlands using remote sensing data
The results obtained from our study demonstrate the successful modeling of small inland wetland occurrence in South-Kivu, eastern DRC, by combining optical and SAR indices using machine learning (ML) algorithms.The  produced wetland map was reclassified into two classes, 'wetland' and 'non-wetland' , and converted into shapefiles.Among the four ML models tested, the Random Forest (RF) model exhibited high accuracy, with an AUC of 0.97 and a TSS statistic of 0.84, indicating strong discrimination between wetlands and non-wetlands areas.These findings align with the results obtained by Garba et al. 32 in Nigeria, where a similar approach was employed, and the study area resembles our environment.Eleven out of the 27 indices were used for the RF model, mainly vegetation indices (WET, MDNWI, GOSSAN, and EVI), SAR ratio (η, ρ), and topographic (TWI, slope, land aspect, curvature, and elevation).MNDWI seems to be a good vegetation index for wetland mapping 32 ; it contributed up to 21%, 23%, 20%, and 18% for ANN, BRT, MaxEnt, and RF, respectively (Table 1).MNDWI is an index currently used for the enhancement of open water features; the index also diminishes built-up area features that are often correlated with open water in other indices (e.g., when using the G and SWIR bands: pixel values from the green and short-wave infrared band respectively) 67 .In South-Africa, Slagter et al. 68 used the same approach in combining both Sentinel 1 and 2 and RF and found that 4 explanatory variables (VV, VH; VV/VH, NDVI, and MNDWI) can be used and advised for wetland mapping.However, the problem persisted consistently in wetland areas characterized by the presence of trees.To address the challenge of mapping wetlands in densely forested areas (such as the Shabunda and Mwenga territories), we anticipated that radar sensors operating in high-wavelength L-or P-band with the HH polarization mode would yield more accurate results, as suggested by Slagter et al. 68 .The HH mode enables the observation of double-bounce scattering during floods, essential for mapping wetlands in highly vegetated areas.However, it was anticipated that C-band sensors operating in the VV polarization mode would have limited capabilities for mapping highly vegetated wetlands 68 .
The observation of double-bounce scattering during floods is required for mapping wetlands in densely forested areas (such as Shabunda and Mwenga territories in our case), and this necessitates a certain degree of vegetation penetration of a radar sensor.This is accomplished using high-wavelength L-or P-band sensors with HH mode in highly vegetated wetlands like mangrove or swamp forests.Varied results have also been obtained with C-band and VV mode 68 .It was anticipated in our study that the C-band sensors, operating in η: VV/VH and γ: HH/VH modes, would have only fair capabilities to map highly vegetated wetlands.In our case, η and ρ ratio contributed 10-12% for RF and 8-10% for ANN.We also anticipated and confirmed that higher-resolution inland wetland mapping with Sentinel and ALOS-PALSAR would capture a smaller wetland size than previously documented in regional datasets 38 .Indeed, integrating the SAR improves the accuracy of models.This can be attributed to the improvement in the accuracy of identifying small wetlands vegetation structures and soil water content captured by SAR imagery.Our findings and from Garba et al. 32 allow us to conclude that the combined use of optical and SAR indices resulted in greater accuracy for small inland wetlands and all wetland classes in general than the use of each one in isolation; such a combination can be advised to produce refined wetland maps.Our final map (Fig. 10) illustrates the spatial distribution of small inland wetlands in South-Kivu province with pixels of 10 m in size.Overall, the classification result shows high accuracy when the RF model is used.This study's post-classification analysis and digitalization revealed that wetlands account for ~ 13.5% (898,690 ha) of the entire province.These wetlands exhibit a wide range of sizes, varying from a few dozen to thousands of hectares.Most of the wetlands are located in lower altitude territories rather than highlands.However, many small wetland fragments, often-single pixels, were observed, particularly in high-altitude areas.These fragments were reclassified using post-classification methodologies outlined in the methodology section.Regarding the geographic characteristics, it is noteworthy that peats, swamps, and ponds may have developed at very high altitudes.Nzabandora and Roche 69 suggest that in the territories of Kabare, Kalehe, and Walungu, as well as in extremely high-altitude regions (> 2700 m) in enclosed valleys, these ecosystems can easly form.In these areas, the average temperature can drop below 10 °C at the top of the Kivu's dorsal, while sufficient light, storms, and rainfall maintain an annual precipitation rate of 1600-1700 mm.The presence of significant fault lines in the Lake Kivu sector of the Congo strongly influences the eastern face of the region.The upper courses of these rivers, situated between 2200 and 2300 m.a.s.l, exhibit a senile appearance and give rise to extensive stretches of wetlands, some of which have developed into peatlands.Within these wetlands, native palustrine vegetation, primarily composed of Cyperus denudatus and Cyperus latifolius, dominates the landscape.These dense vegetation cover in certain areas might explain why traditional vegetation indices, such as the NDVI, SAVI, etc., did not significantly impact the wetland mapping process.Nevertheless, gaining a comprehensive understanding of wetland fragmentation in eastern DRC and its impacts on biodiversity, ecosystem services provided to the community, and the role of both larger and smaller inland wetlands in the regional landscape will require further research.The finding from this dataset serves as a valuable starting point for future modeling efforts aimed at enhancing our understanding of these effects.Figures 4, 11 and 13 provide a closer view of some wetland complexes on the wetland map, highlighting their significance despite variations in size and surface area.
The comparison of the produced wetlands map with existing maps in the vicinity revealed interesting insights.It was found that the extent of wetlands in South-Kivu province, as depicted in our study, was more significant compared to some previous studies, such as the one conducted by Kulimushi et al. 46 and the map available at the Center for International Forestry Research (CIFOR).The CIFOR map was generated at a spatial resolution of 231 m, employing transparent rules on hydrological wetness, satellite-derived soil wetness phenology, and geomorphology.However, the wetlands' extent was underestimated due to its global scale.In contrast, our proposed methodology integrated three biophysical indices that capture essential wetlands characteristics, namely, (i) long-term water supply exceeding atmospheric water demand; (ii) annually and seasonally waterlogged soils; and (iii) a favorable geomorphological position for water provision and retention Gumbricht 70 .
Nonetheless, there is still room for improvement in our methodology.For instance, the inclusion of SAR data captures seasonal variability and enhances the accuracy of wetlands mapping.However, using SAR data is often limited due to their high cost and complex processing requirement 67,71,72 .Nevertheless, recent initiatives such as Sentinel-1 and -2 and ALOSPALSAR have made radar data more accessible for wetland studies.Delineating wetlands through in-situ examination of hydric soil characteristics is time-consuming and expensive, as previously mentioned in our introduction and Lidzhegu et al. 67 .As an alternative, optical remotely sensed www.nature.com/scientificreports/images face challenges such as cloud cover and spectral confusion among different land cover categories.To address these issues, incorporating SAR data can penetrate through clouds and provide valuable information.However, providing these data is sometimes possible thanks to Sentinel 1 and 2 and ALOSPALSAR.However, the computational complexities of processing SAR data should be considered 67,71,72 .Regarding the bands used in radar systems, different frequencies are employed for wetland mapping.These include P-band: ~ 69.0 cm (BIO-MASS), L-band: ~ 23.5 cm (ALOS-2 PALSAR-2, SAOCOM-1, NISAR-L), S-band: ~ 9.4 cm (NovaSAR, NISAR-S), C-band: ~ 5.6 cm (Sentinel-1, Radarsat-2, RCM), X-band: ~ 3.1 cm (TerraSAR-X, TanDEM-X, COSMO-SkyMed).This study used C-and -L bands for ALOSPALSAR and Sentinel since they are freely available compared to other SAR data 72 .While our approach did not integrate geomorphological and pedologic variables, we found that topographic indices played a significant role in wetland distribution.Variables such as slope, TWI, Curvature, and other indices like TRI and SPI contributed substantially to the wetland mapping.The compound topographic index (CTI), also referred to as the topographic wetness index (TWI), is particularly relevant as it quantifies steadystate moisture and exhibits strong correlations with various soil characteristics 67,72 .Additionally, Lidzhegu et al. 65 and Ludwig et al. 73 noted that although numerous spectral indices for water and wet soil detection exist, wetlands can still be easily confused with other land cover types, such as forest and shadows, due to similar spectral profiles.Therefore, our approach incorporated additional auxiliary data in the classification process to mitigate such errors.Previous studies have also highlighted the influence of topographic indices on wetlands mapping.Indeed, Hansen et al. 74 , Rapinel et al. 11 , Guasselli et al. 75 ; Berhanu et al. 56 utilized digital terrain models (DTMs) as explanatory variables to identify estuarine-fringe wetlands.Three topographic variables derived from DTMs were commonly used: (a) the multiscale topographic position index (TPI), (b) the vertical distance to the channel network (VDCN), and (c) the topographic wetness index (TWI), TWI, which characterizes potential soil wetness based on contributing area and local slope, typically ranges from 0 to 30, with higher values indicating a higher probability of wet soil 67 .TWI has proven to be suitable for characterizing riverine wetlands.We integrated two indices, mainly TWI and TPI, which significantly contributed to the wetland mapping process.
The TWI also helps predict soil characteristics, including horizon depth, silt percentage, organic matter content, and phosphorus.The index is used to define biological processes such as yearly net primary production, vegetation patterns, and forest site quality, as well as to explore spatial scale effects on hydrological processes and identify hydrological flow pathways for geochemical modeling 56,75,76 .However, according to Ludwig et al. 73 , even though a variety of spectral indices for water and wet soils detection (hereafter referred to as wetness) is available, wetlands are still easily confused with other upland land cover types such as forests and shadows, since they share similar spectral profiles.In our approach, we have included additional auxiliary data in the classification to minimize these errors significantly.

Accuracy of ML models
According to the literature review, the four models tested for wetland mapping are commonly used in central and East Africa studies.Among these models, Random Forest 'RF' demonstrated higher accuracy than the other three.Only the ANN model comes close to the accuracy of RF.This suggests that RF is the most accurate model for mapping small inland wetlands, and its capabilities in wetland mapping and monitoring have been consistently demonstrated.This conclusion aligns with the findings of Garba et al. 32 in Nigeria, Barbosa and Maillard 77 in Brazil, and Slagter et al. 68 in South-Africa, who also identified RF as a highly accurate model for wetland mapping.Slagter et al. 68 further discussed the integration of SAR data and found no significant accuracy differences between Sentinel-1 and -2 for mapping surface water dynamics.Our results are similar to Slagter et al. 68 , as we also observed high accuracy when combining Sentinel-1 and -2 data.Whyte et al. 33 conducted a study in South Africa using both Sentinel-1 and -2, as well as SVM and RF models, with RF demonstrating higher accuracy (83.3% OA, Kappa = 0.72) compared to SVM (79.8% OA, Kappa = 0.68).These accuracies were lower than those achieved with optical data alone but increased when optical and radar data were combined.
Despite using the same RF model, the difference in accuracy between our study and Whyte et al. 33 could be attributed to several factors.One possibility is the integration of new vegetation indices in our study or the specificity characteristics of our region.It is important to note that complex data patterns can be unique to specific geographies, and a model trained in one geographic landscape may not perform equally well in different geography 33 .In our study, the selection of parameters for the four classifiers (ANN, BRT, Maxent, and RF) allowed for a fairer comparative analysis rather than relying on specific classifier evaluations.This technique has been successfully implemented in other LULC investigations [78][79][80] .Across all evaluation metrics, RF consistently outperformed ANN, BRT, and Maxent, as demonstrated by the AUC, Kappa, correlation, and TSS values.The statistical tests conducted, such as the DeLong Test, F welch, and Brown-Forsythe tests, further confirmed the significant superiority of RF over the other models (as shown in Fig. 4).
The observed differences in the lowest user accuracies between RF and other models can be attributed to the processing steps of each model.At the same time, both RF and BRT combined decision trees, and BRT started the combination process earlier.The RF help to reduce the variation observed in decision trees: by (i) employing various training samples, (ii) defining sub-ensembles with random characteristics, and (iii) building and combining shallow trees (slightly deep) 59 .
Many potential reasons for the different results of the four models can be mentioned.In fact, in each of these classifiers, there are factors that can affect the classification accuracy.Among these, the image segmentation step, the training sample and feature selection, and the parameters tuning set as advised by Mahdavi et al. 37 .While McNairn et al. 81 and Adam et al. 33 mentioned that classification accuracy is not the only thing to be considered on a classifier, the operational monitoring purposes, use friendly, the objective of the study and the type of the study area (here small wetlands), etc. have also to be considered.Nevertheless, these classifiers are among the most famous classification algorithms used for wetland mapping 37 .BRT is a supervised classifier belonging to classification and regression trees (CART), here input data are divided into mutually exclusive groups based on their attributes which is different from RF, a Bayesian statistic assuming that feature vectors of each class are normally distributed.While considered as a CART, RF is considered as an extension of Decision tree (DT).Nevertheless, each of these models has its own strengths and weakness.They differ in term of structure, composition and learning process.The choice should be based on the specific characteristics of the available data, the goals of analysis and resources available.Based on the model performance, the specific of the small wetland context, and the quality of our input data, RF was selected as the best model.Supplementary data 14 presents the strengths and weakness identified for each model used.In summary, RF is an ensemble learning algorithm suitable for various tasks, while Maxent is specifically designed for species distribution modeling, and ANN is a versatile algorithm capable of handling complex patterns in data.
This study presents an affordable and practical technique for accurately delineating small inland wetlands using freely available data at a reasonable spatial resolution.Although Landsat is commonly used for wetland mapping in eastern Africa 82 , we opted for Sentinel-1 and -2 and ALOSPALSAR, despite their lower resolution compared to purchased satellite images from providers such as WorldView, Pleiades, GeoEye, etc. which remain very expensive and almost impossible to obtain at the scale of the entire province.These satellite data sources allowed us to achieve satisfactory accuracy in mapping wetlands.
One advantage of our methodology is that we utilized RStudio packages and freely available scripts, eliminating the need for expensive software.This makes the technique accessible and cost-effective for researchers and practitioners involved in wetland mapping.However, it is essential to note that the RF and ANN required significantly more computation time than other classifiers.In our case, using a desktop computer with Intel(R) core (TM) i7-11800 h CPU with 3.5 GHz processor and 32 Go of RAM, the RF and ANN computations took approximately 22-24 h to complete.This longer processing time could challenge more extensive and longterm studies.More powerful hardware with faster processors and greater RAM capacity could be employed to overcome this issue.Furthermore, the execution time could be reduced by optimizing the algorithm implementation or utilizing more powerful hardware.Including additional variables for monitoring purposes, such as multi-temporal data spanning several years, may also affect processing time.Another potential solution is to explore using the Google Earth Engine (GEE), a cloud-based platform offering extensive geospatial processing capabilities.GEE has the advantage of scalability and efficiency, enabling parallel computing and the integration of various data sources.By leveraging GEE, computational challenges associated with large-scale wetland studies can be addressed, opening opportunities for further development and refinement of the methodology.Despite the computational considerations, our study demonstrates the feasibility of cost-effective wetland mapping using freely available data.As technology continues to advance and more powerful computing resources become available, the efficiency of the process can be improved, facilitating larger-scale studies and supporting ongoing wetland monitoring efforts.

Study limitation
While this study has contributed significantly to our understanding of the spatial distribution of small inland wetlands in the South-Kivu province, it does have certain limitations that should be acknowledged.One limitation is related to the accessibility of certain territories within the province.Due to the challenges of accessing some areas, fewer training sites were established than in other territories.This discrepancy in training samples could result in variations in precision and accuracy across different territories.Territories closer to urban centers and more accessible, such as Kabare, Walungu, Uvira, and Kalehe, had more training points than Mwenga, Shabunda, and Fizi.This uneven distribution of training samples could introduce bias and affect the accuracy of wetland mapping in different regions.Future research should address this limitation by ensuring a more balanced representation of training samples across all territories.
Another limitation of the methodology used in this study is grouping all wetlands into a single class.While the images in Fig. 4 demonstrate the diversity of wetland types, ranging from lakes and river shore wetlands to peatlands, marshlands, bogs, and swamps, the available training samples for these specific wetland types were limited.For instance, only a small number of peatland samples were identified during fieldwork (~ 12 only).As a result, all these wetland classes were merged to reduce classification errors.Therefore, future research should aim to map and differentiate between these various wetland types to improve wetland maps' accuracy and representational quality.
Additionally, the study lacks lithological, geological, and pedological variables.The available data on these aspects have a relatively low spatial resolution, ranging from 5 to 15 km 83 .Despite the potential relevance of these variables in delineating wetland types, they were not integrated into the analysis due to their limited resolution.However, indices such as TRI and TWI, which are associated with these elements, have been included.It is important to mention that the available soil data have resolutions ranging from 500 m to 5 km.However, we aimed to restrict the resolution of our results to 10 m.Therefore, this decision also justified excluding these data from the integration process.
Furthermore, this study primarily focused on mapping and delineating small inland wetlands without providing a comprehensive characterization of these wetlands, their ecosystem services, and constraints in utilization.To fully understand the extent and functioning of different wetland types and their fragmentation, future research should consider comprehensive assessments incorporating detailed information on wetland characteristics, ecosystem services provided, and the impacts of human disturbances.Such comprehensive wetland mapping and knowledge of their fragmentation patterns are crucial for economic assessments and decision-making by regional and international agencies.Also, while the study used images from two different seasons to assess the seasonality of inundated areas, there are limitations in capturing the peak of inundations 33 .The periodicity of satellite imagery may not align with the specific timing of peak inundation events.Additionally, the strong interannual variability of African river and water flow regimes poses challenges in accurately detecting inundation extents using publicly available satellite data.Moreover, cloud cover can sometimes hinder the visibility of wetland areas during satellite image acquisition.Addressing these challenges would require the development of robust cloud masking models and potentially exploring alternative data sources or techniques, such as synthetic aperture radar (SAR), to overcome these limitations.
In summary, while this study has advanced our understanding of small inland wetlands in the South-Kivu province, it is essential to acknowledge the limitations related to the distribution of training samples, the grouping of wetlands into a single class, the exclusion of certain environmental variables, and the challenges in capturing peak inundation events and addressing cloud cover [84][85][86] .Future research should strive to address these limitations and incorporate a more comprehensive approach to wetland mapping and characterization.

Conclusion
Based on the findings of this study, we can confidently conclude that the mapping of small inland wetlands in the South-Kivu province can be carried out with remarkable precision using a combination of topographic and vegetation indices within a Random Forest (RF) model.By incorporating Synthetic Aperture Radar (SAR) data, we were able to enhance the accuracy and effectively capture the seasonal variations of these wetland areas.Our proposed methodology, which involved carefully selecting a subset of variables and considering new indices, yielded an impressive accuracy rate of approximately 72%.Notably, variables such as η and ρ backscattering ratio, MDNWI, TWI, slope, and elevation played a significant role in achieving these results.Our analysis estimated that ~ 13.5% (equivalent to 898,690 ha) of the South-Kivu province is covered by small inland wetlands.These wetlands exhibit a wide range in size, spanning from a few acres to vast expanses of thousands of hectares.It is important to acknowledge that due to data limitations, we merged different types of wetlands into a single class to avoid introducing bias.Nonetheless, this approach creates exciting opportunities for future research endeavors to delve into the characterization and classification of the diverse wetland types identified during our fieldwork.The significance of wetlands in vital ecological processes such as the water cycle, greenhouse gas exchange, carbon dynamics, and the support of biodiversity cannot be overstated.As such, our study serves as a stepping-stone for further investigations aimed at comprehending the functionality, ecosystem services, and potential risks associated with these wetlands.
Future research should focus on assessing the specific ecosystem services provided by these wetlands, quantifying their contribution to carbon storage, and evaluating their role in supporting and preserving biodiversity.Our findings underscore the importance of understanding and conserving wetland ecosystems, they provide valuable insights for informed decision-making regarding wetland conservation, management, and sustainable land use practices in the South-Kivu province.By understanding these unique ecosystems comprehensively, we can effectively protect their invaluable services, mitigate risks such as floods and pollutants, and harness their potential for climate change mitigation.

Figure 2 .
Figure 2. Diversity of wetlands in South-Kivu, Eastern DR Congo (a) very high altitude lake Lubwe in Itombwe, (b) peatland in Lugana, (c) coastal of the Ruzizi river, and (d) flooded plain for rice production in Ruzizi plain.map created using ArcGIS 10.7 Esri-TM: http:// www.esri.com.

Figure 3 .
Figure 3. Flowchart of the methods used to map wetlands in South-Kivu by using optical, SAR and field data and four machine-learning (ML) algorithms.Topographic, vegetation, hydrologic indices were used as variables subsequently calibrated and validated using field samples.

Figure 4 .
Figure 4. Model performances using AUC, Kappa and TSS, correlation coefficients.The four models were compared with the F-test of Welch and DeLong test.For each model 10 experiments were executed and helped for statistical comparison.

Figure 5 .
Figure 5. Images obtained from fieldwork when identifying and delineation wetlands in South-Kivu province, here: (a) swamp in Kabare, (b) papyrus wetland along the Lake Kivu in Kalehe and Kabare, (c) marshland in Kabare, (d) and (e) inundated valley in Walungu, (f) peatland in the Nakananda wetland in Walungu, rice valley in Ruzizi plain.

Figure 6 .Figure 7 .
Figure 6.Prediction accuracy of the four models with and without integration of SAR data (here only the AUC and TSS were integrated as accuracy indices).

Figure 8 .
Figure 8. Proportion of wetlands and non-wetlands in the South-Kivu province (a) and in the eight territories of South-Kivu province (b), (Here we consider only results from RF when we combined optical and SAR data.(a) Was obtained after logarithm "log2" conversion of x-axis values for better visualization).

Figure 9 .
Figure 9. Surface (ha) and perimeter (in km) of the two classes (wetland and non-wetland) after clipped for each territory scale (the eight territories were split here in large and small territory for data representation and visualization.

Figure 10 .
Figure 10.Map of wetland produced in the South-Kivu province, eastern DR Congo using the four models (a), and permanently and seasonally flooded wetland map obtained after the integration of SAR data from wet and dry season in RF model (b).(Maps created using ArcGIS 10.7 Esri-TM: http:// www.esri.com).

Figure 12 .
Figure 12.Wetland area in the eight territories of the South-Kivu province, eastern DRC (values were transformed in log2 for better visualization.

Table 1 .
Variable importance for both SAR and Optical indices delivered from Sentinel and ALOS-PALSAR.